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Abstract — A novel subcarrier-pair based opportunistic DF pro- 
tocol is proposed for cooperative downlink OFDMA transmission 
aided by a decode-and-forward (DF) relay. Specifically, user 
message bits are transmitted in two consecutive equal-duration 
time slots. A subcarrier in the first slot can be paired with a 
subcarrier in the second slot for the DF relay-aided transmission 
to a user. In particular, the source and the relay can transmit 
simultaneously to implement beamforming at the subcarrier in 
the second slot for the relay-aided transmission. Each unpaired 
subcarrier in either the first or second slot is used by the source 
for direct transmission to a user without the relay's assistance. 
The sum rate maximized resource allocation (RA) problem is 
addressed for this protocol under a total power constraint. It is 
shown that the novel protocol leads to a maximum sum rate 
greater than or equal to that for a benchmark one, which 
does not allow the source to implement beamforming at the 
subcarrier in the second slot for the relay-aided transmission. 
Then, a polynomial-complexity RA algorithm is developed to 
find an (at least approximately) optimum resource allocation 
(i.e., source/relay power, subcarrier pairing and assignment to 
users) for either the proposed or benchmark protocol. Numerical 
experiments illustrate that the novel protocol can lead to a much 
greater sum rate than the benchmark oneQ 

I. Introduction 

The incorporation of subcarrier-pair based decode-and- 
forward (DF) relaying into orthogonal frequency division mod- 
ulation (OFDM) or multiple-access (OFDMA) transmission 
and associated resource allocation (RA) were studied in |IT]- 
m when the source-to-destination (S-D) link exists. In |[T]-|l3l, 
an "always-relaying" DF protocol was used, i.e., a subcarrier 
in the first time slot is always paired with a subcarrier in 
the second slot for the relay-aided transmission. To better 
exploit the frequency-selective fading, we have proposed an 
opportunistic DF relaying protocol (sometimes termed as 
selection relaying) in ||4|-||6l, i.e, a subcarrier in the first time 
slot can either be paired with a subcarrier in the second slot 
for the relay-aided transmission, or used for the S-D direct 
transmission without the relay's assistance. It is very important 
to note that when some subcamers in the first slot are used 
for the direct transmission, some subcarriers in the second slot 
will not be used, which wastes spectrum resource. 

To address the above issue, we have proposed an improved 
DF protocol in |7j. This protocol is the same as those 

' Research supported by The Program for Professor of Special Appointment 
(Eastern Scholar) at Shanghai Institutions of Higher Leaming. It is also sup- 
ported by the TAP project BESTCOM, the ARC SCOOP and the NEWCOM# 
and NSF China #61271213. 



considered in ID except that the source can also make direct 
S-D transmission at every unpaired subcarrier in the second 
slot. This protocol and its RA were later intensively studied, 
e.g., in II7II- III2I . Note that the improved protocol does not 
really improve the way how DF relaying is implemented over 
a subcarrier pair, but rather let the source utilize the unpaired 
subcarriers in the second slot for direct transmission to avoid 
the waste of spectrum resource. In ini-ini, the subcarrier 
pairing and power allocation are jointly optimized for point- 
to-point OFDM transmission. As for OFDMA systems, RA 
problems considering the joint optimization of power alloca- 
tion and subcarrier assignment to users are addressed in It7|- 
ifTOl . In these works, a priori and CSI-independent subcarrier 
pairing is considered, i.e., a subcarrier in the first slot is always 
paired with the same subcarrier in the second slot if the relay- 
aided mode is used. It is a complicated RA problem to jointly 
optimize subcarrier pairing, power allocation and subcarrier 
assignment to users. 

In this paper, we consider downlink OFDMA transmission 
from a source to multiple users aided by a DF relay. Compared 
with the existing works, this paper makes the following 
contributions. First, a novel subcarrier-pair based opportunistic 
DF relaying protocol is proposed. A benchmark protocol using 
the improved protocol as in Q is also considered. How- 
ever, the proposed protocol uses further improved relay-aided 
transmission, which allows the source and relay to transmit 
simultaneously to implement beamforming at the subcarrier 
in the second slot. Note that the proposed protocol truly 
improves the implementation of DF relaying over a subcarrier 
pair with transmit beamforming, which is not the case for 
the benchmark protocol. Second, the sum rate maximized 
RA problem is addressed for both the novel and benchmark 
protocols, under a total power constraint for the whole system. 
It is shown that the novel protocol leads to a maximum sum 
rate greater than or equal to that for the benchmark one. An RA 
algorithm is developed for each protocol to find the globally 
optimum source/relay power allocation and subcarrier pairing 
to maximize the sum rate of all users. 

The rest of this paper is organized as follows. In the next 
section, the system and transmission protocols are described. 
In Section IIII-AI we will focus on computing the maximum 
rate and optimum power allocation for a subcarrier pair using 
the relay-aided transmission for both protocols. Using these 
results, an RA algorithm will be developed in Section |III1 and 



numerical experiments are shown in Section |IV] Finally, some 
conclusions are drawn. 

Notations: A letter in bold, e.g. x, represents a set. C{x) = 
\ log2(l + x) and [x\^ = max{x, 0}. 

II. Protocols and WSR maximization problem 
A. The transmission system and protocols 

Consider the downlink OFDMA transmission from a source 
to U users collected in the set U = {u\u = 1, • • ■ ^U} aided 
by a DP relay. The source and the relay can simultaneously 
emit OFDM symbols using K subcarriers and with sufficiently 
long cyclic prefix to eliminate inter-symbol interference. User 
message bits are transmitted in two consecutive equal-duration 
time slots, during which all channels are assumed to keep 
unchanged. During the first slot, only the source broadcasts 
N OFDM symbols. Both the relay and all users receive these 
symbols. After proper processing explained later, the source 
and relay simultaneously broadcast N OFDM symbols, and 
the users receive them during the second slot. Due to the 
OFDMA, each subcarrier is dedicated to transmitting a single 
user's message exclusively. A subcarrier in the first slot can 
be paired with a subcarrier in the second slot for the relay- 
aided mode transmission to a user Each unpaired subcarrier 
in either the first or second slot is used by the source for the 
direct mode transmission to a user 

To simplify description, we use subcarriers k and I to denote 
the fcth and Zth subcarriers used during the first and second 
slots, respectively {k,l = 1, • " ' iK)- We define the source 
transmission powers for subcarrier k in the first slot and 
subcarrier I in the second slot as Ps.fc.i and Ps,i.2, respectively. 
The relay transmission power for subcarrier I is Pr,i.2- The 
complex amplitude gains at subcarrier k for the source-to- 
relay, source-to-u and relay-to-u channels are /isi.fc, ^su.fc and 
hvu,k, respectively. The two transmission modes for the novel 
protocol are elaborated as follows: 

1} The relay-aided transmission mode: Suppose subcarrier 
k is paired with subcarrier I for the relay-aided mode transmis- 
sion to user u. A block of message bits are first encoded into 
a code word of complex symbols {9{n)\n = 1, • • • , N} with 
E{\0{n)\'^) = 1, V n. In the first slot, the source broadcasts the 
codeword over subcarrier k as illustrated in Figure [T]a. At the 
relay and user u, the nth baseband signals received through 
subcarrier k are 



and 



yu,k,i{n) \/ Ps,k,ihsu,k6{n) + 2:„,fc,i(ri), n = 1, 



(2) 



respectively, where Zi-.fc(n) and Zu.k,i(n) are both additive 
white Gaussian noise (AWGN) with power a^. The signal- 
to-noise ratio (SNR) at the relay is Ps,k.iGsr,k where Gsr,fe = 

1^ |2 

' ""'J' . At the end of the first time slot, the relay decodes 
the message bits from {y-r^k{n)\n = 1, • • ■ ,N} and then 




lubcarrier k in the fiist slot 



lubfiii-iier I in the second slot 



(a) 



(b) 



Fig. 1. The relay-aided transmission mode over the subcanier pair {k, I) to 
user u. 



reencodes those bits into the same codeword as the source 
did. 

In the second time slot, the source and relay broadcast 
the codewords {6'(n)e~^^''=">' |V n} and {6'(n)e"-'^''-^" ' |V n} 
through subcarrier /, respectively, where Zhsu.i and Zhyu,i 
represent the phase of hsu.i and hjuj, respectively. This means 
that the source and relay implement transmit beamforming to 
emit the codeword through subcarrier I as illustrated in Figure 
[T]b. Note that the source and relay need to know the phase 
of hsu.i and hmj, respectively. At user u, the nth baseband 
signal received through subcarrier I is 



y«,i,2('^) = {\/ Ps,l,2\^su,l \ + ^/P^\hl■u,l\)0{n) + Zu,l.2{n), 

(3) 

where Zu.i,2{'n) is the AWGN with power cr^. 

Finally, user u decodes the message bits from all signals 
received during the two slots. These signals can be grouped 
into N vectors, the nth of which is 



y(") = 



where z(n) 



Ps,k,lhsu,k 
\/PsJa\hsuJ 



(4) 



6'(n) +z(n), 



[zu.,kA{n),Zu,i.2{'>T)] ■ Note that the transmis- 
sion of the codeword in effect makes N uses of a discrete 
memoryless single-input-two-output channel specified by (|4|, 
with the nth input and output being 6'(n) and y(n), respec- 
tively. To achieve the maximum reliable transmission rate, 
maximum ratio combining should be used ||141 . It can readily 
be derived that the SNR for after this combining is 



where G. 



'yklu{Ps,k,l, Psd,2, Pi,l.2) = Gsu,kPs,k,l + 
{\/ GsujPs,l.2 + \/G 

- and G,„, = 



(5) 



To ensure both 



the relay and user u can reliably decode the message bits, 
the maximum number of message bits that can be transmit- 
ted is 2NCiGsrMPsMs) and 2NC{jklu{PsMS, Ps,l,2, PrJ,2)), 

respectively. This means that the maximum transmission rate 
over the subcarrier pair (fc, I) in the relay-aided mode to user 

u is equal to C{mm{Gsi,kPs,k,i,^kiu{Ps.k.i, Ps,L2, P1J.2)}) 
bits/OFDM-symbol (bpos). 

2) The direct transmission mode: Suppose subcarrier k 
(respectively, subcarrier /) is unpaired with any subcarrier in 
the second (respectively, first) slot, and is used for direct mode 
transmission to user u. The source first encodes message bits 



into a codeword of N symbols, which are then broadcast 
through subcarrier k (respectively, subcarrier I. In such a case, 
the relay keeps silent at subcarrier I in the second slot, i.e. 
^r./,2 = 0.). User u decodes the message bits from the signals 
received through subcarrier k (respectively, subcarrier I). The 
maximum rate through subcarrier k (respectively, subcarrier I) 
in the direct transmission mode is C{Ps,k.iGsu,k) (respectively, 
C{Psd,2Gsu,i)) bpos. 

A benchmark protocol is also considered. This protocol 
is the same as the novel protocol except for the relay-aided 
transmission mode. Specifically, the relay-aided mode is the 
same as that widely studied in the literature, i.e., the source 
does not transmit at subcarrier I during the second slot, if 
subcarriers k and I are paired for the relay-aided transmission 
to user u. In such a case, the maximum rate for the relay- 
aided transmission over that subcarrier pair to user u is equal 
to C( 

B. The sum rate maximization problem 

We assume there exists a central controller which knows 
precisely the CSI {Gsr:.k,Gsu.k,Gi-u.k\^ k}. Before the data 
transmission, the controller needs to find the optimum sub- 
carrier assignment and power allocation to maximize the sum 
rate of all users for the adopted transmission protocol (which 
can be either the novel or benchmark protocol), when the total 
power consumption is not higher than a prescribed value Pt- 
Then, the controller can inform the source and the relay the 
optimum RA to be adopted for data transmission. 

It can be shown that the novel protocol leads to a maximum 
sum rate greater than or equal to that for the benchmark 
one. This can be proven as follows. Note that the benchmark 
protocol is a special case of the novel one, since it is equivalent 
to the novel one constrained with Ps.1,2 = if subcarrier I is 
paired with a subcarrier for the relay-aided transmission. Sup- 
pose the optimum subcarrier assignment and power allocation 
have been found for the benchmark protocol. By using the 
novel protocol with the same subcarrier assignment and power 
allocation, the same sum rate can be achieved. Obviously, the 
maximum sum rate for the novel protocol is greater than or 
equal to that sum rate, namely the maximum sum rate for the 
benchmark protocol. 

III. RA Algorithm design 

A. Rate maximization for the relay-aided mode over a sub- 
carrier pair 

Assume subcarriers k and I are paired for the relay-aided 
mode transmission to user u, and a sum power P is used 
for this pair. For the benchmark protocol, it can be shown 
by using an intuitive method similar as the one introduced in 
Appendix of [91 that the maximum rate associated with the 
above optimum solution is equal to C{G\i.^P) with 



G 



if min{Gsr,fe,G,.„,i} > Gsu,k, 
min{Gsr,fc,Gs„,fc} if min{Gsr,fc, G™,;} < Gs„,fc, 

(6) 



and the optimum Ps,fc,i is 



P. 



k.l 



A„ fc+G,„ min{Gsr,fe,G™,;} > Gsu^fc, 

P if mm{Gsr,k,G,ud} < G^u^k- 



To facilitate the derivation for the proposed protocol, define 

A„.fe = Gsr,fc - Gsu.k and Gu,i = Ggu,; +G,.„,;. To maximize 
the rate, the optimum Ps.k.i^ Ps,i.2 and Pr,i.2 are the optimum 
solution for 

max mm{Gsr,kPs,ks,lkiu{Ps,k,i, Ps,L2, Pi,i,2)} 

S.t. P.^kS + Ps,U2 + PrJ,2 = P, (7) 
PsMS >0,Ps,i,2 >0,PrJ,2 >0. 

Using Cauchy-Schwartz inequality and an intuitive method 
similar as that introduced in Appendix of f9], it can be shown 
that the optimum Pg.k.i^ Ps,i.2 and Pr.i,2 for Q are 



Ps 



fe,i 



P 



^P if min{Gsi.,fc,Guj} > Gsu^fc, 
if min{Gsi.,fc,Gu^;} < Gsu.k, 



P. 



1,2 



G^'i (Au k+Gu i) ^ niiniCsr^fci Gu^;} > Gsu,ki 

if 



and Pr,i.2 = P — Ps,k,i — Ps.i,2 (please see fTSl for more de- 
tails). The maximum rate associated with the above optimum 
solution is equal to C{G^i^P) with 



Gn 



Gsr,kGu,l 



if min{Gsr,fc, G„,/} > G 
min{Gsr,fc,Gs„,fc} if min{Gsr,fc, G„,/} < Gs 



(8) 



B. Formulation of the WSR maximization problem 

For the adopted protocol (which can be either the proposed 
or benchmark protocol), we first define 



Gkiu 



Gkiu if fl^*^ proposed protocol is adopted, 
Gtiu fi^*^ benchmark protocol is adopted. 



For any possible subcarrier assignment used by the adopted 
protocol, suppose m subcarrier pairs are assigned to the relay- 
aided transmission, the unpaired subcarriers in the two slots 
can always be one-to-one associated with each other to form 
K — m virtual subcarrier pairs for the direct transmission. 
Based on this observation, we define: 

• tkiu G {0,1} and Pkiu > 0, V k,l,u. tkiu = 1 indicates 
that subcarrier k is paired with subcarrier / for the relay- 
aided transmission to user u. When tkiu = 1. Pkiu is 
used as the total power for the subcarrier pair (fc, /). 

• tkiab G {0,1}, akiab > and fikiab > 0, V k,l,u. 
tkiab ~ 1 indicates that subcarrier k is assigned in the 
direct transmission mode to user a during the first slot, 
and so is subcarrier / to user b during the second slot. 
When tkiab = 1, PsMS and Psj,2 take the value of akiab 
and Pkiab, respectively. 



Let us collect all indicator and power variables in the sets I 
and P, respectively, and define S = {I, P}. Every feasible RA 
scheme can be described by an S satisfying simultaneously 



tkiu,tkiab e {0, 1}, V fc, l,u, a, h, 

tkiu + y] tkiab =i,v k, 



(9) 
(10) 



J2{J2^kiu + Y.tkiab\^hyi, (11) 

k \ u a,b j 

{tkluPklu + tklab{oLklab + Pklab)) < ^t, (12) 

Pkiu > 0, akiab > 0, 13kiab > 0, V fc, M, a, 6. (13) 

Given a feasible S, the maximum sum rate for the adopted 
protocol is 

/(S)= J2 {tkiuCiGkiuPkiu)+ (14) 

k,l,u,a,b 

tklab{C{Gsa,kCy.klab) + C{Gsb,lPklab)) , 

and the sum rate maximization problem is to solve 

max /(S) (15) 
s.t. ©-O 



for a globally optimum S. Obviously, (flSl l is a nonconvex 
mixed-integer nonlinear program. To find a globally optimum 
S, all indicator variables are first relaxed to be continuous 
within [0, 1], i.e.. Then, we make the COV from P to P = 

{Pkiu,akiab, /3feia&|Vfc, /, u, a, b}, where every Pkiu, akiab and 
I3kiab satisfy, respectively, 

Pkiu = tkluPklu, OLklab = tklabCtklab, l^klab = tklabPklab- 

After collecting all variables into X = {I,P}, the RA 
problem can be rewritten as 

max ,g(X) (16) 

s.t. tkiu,tkiab k,l,u,a,h. (17) 

X] {Pklu + auab + Pklab) <Pt, (18) 
k,l,u,a,b 

Pkiu > 0, akiab > O^Pkiab > 0,Vfc, /,it,a, 6, (19) 
where g(X.) represents the maximum sum rate expressed as 
ff(X)= ^ {Htkiu,Pkiu,Gkiu) (20) 

kdjU.a.b 

+ <l>{tklab, Oiklab, Gga.k) + 0(^fciob, l^klab, Gebd)) , 



and 



t-C(Cf) ift>0, 
if t = 0. 



(21) 



Obviously, (fTSI l is a relaxation of ( fTsl l. We will find an 
(at least approximately) optimum solution for (fT6t . and show 
that the S corresponding to this solution is still feasible, and 



hence (at least approximately) optimum for (fTsl l. which can 
be shown in a similar way as reported in lITSl . To this end, 
note that that (f>{t,p, G) is a continuous and concave function 
if t > and x, because it is a perspective function of C{Gp) 
which is concave of p lfT6l . As a result, (/(X) is a concave 
function of X in its feasible domain for ( fTSI l. This means 
that ( fTSI l is a convex optimization problem. Apparently, it also 
satisfies the Slater constraint qualification, therefore its duality 
gap is zero, which justifies the use of the dual method to look 
for the globally optimum of (fTST i, denoted by X* hereafter. 

To use the dual method, i^l is introduced as a Lagrange 
multiplier for the constraint ([18). The Lagrange relaxation 
problem (LRP) for ([16) is 

max X) = .g(X) + M (^^t - P(X)^ (22) 

s.t. dni, (doll - ([ii]i, <0, 

where X) is the Lagrangian of ( fT6] l and P(X) is the 
sum of all Pkiu, otkiab and Pkiab in X. A global optimum 



of (|22| | is denoted by X^. The dual function is defined as 
di^jj) = L(/x,X^). In particular, Pt — -P(Xp) is a subgradient 
of d{ij), i.e., it satisfies V ^',d{fi') > d{fi) + {fi' — n){Pt — 
P(Xp)), and the dual problem is to find the dual optimum 
11* = argmin^>o d{p). 

Since ( fT6l ) has zero duality gap, it satisfies two important 
properties. First, p* > 0. This is because /i* represents the 
sensitivity of the optimum objective value for ( fT6] l with respect 
to Pt, i.e., ^21^ ^* Obviously, g(X*) is strictly 

increasing of Pt, meaning that /i* > 0. Second, /i and X 
are equal to p ~ p* and X^ = X*, if and only if X,, is 
feasible and /i(Pt — P(X^)) = is satisfied according to 
Proposition 5.1.5 in ifTTl . Based on the above properties, the 
fi > and Xp satisfying P(X^ ) = Pt can be found as p* and 
X*. Therefore, the key to developing a duality based algorithm 
consists of two procedures to find fi* and X^, respectively. We 
first introduce the one to find X^ as follows. 

1 ) Finding X^ when ji > Q: The following strategy is used 
to find X^ for ( |22] | when p > 0. First, the optimum P for 



with fixed I is found and denoted by Pj. Define Xj = 
{I, Pi}. Then we find the optimum I to maximizing Xj) 
subject to iT% . ( fTOb and ( fTTI ). Finally, Xj corresponding to this 
optimum I can be taken as X^. 

Suppose I is fixed, we find Pj as follows. Specifically, 
every Pkiu in Pi is equal to when tkiu = 0. When 
tkiu > 0, the optimum Pkiu can be found by using the 
KKT conditions related to Pkiu- In summary, the optimum 
Pkiu can be shown to be Pkiu = tkiuHl^,Gkiu), where 
A(/i, G) — — . In a similar way, the optimum 

ctkiabjind Pkiab can be shown to be auab = tkiabMfJ., Gsa,k) 
and Pkiab = tkiabM.Wb,p,Gsb,i)- respectively. Using these 
formulas, Xi = {I, Pi} can be found. It can readily be shown 
that 

L{p,Xi) = fJ.Pt + y^ {tkluAklu + tklabBuab) (23) 
k.l.u.a.b 



where 



and 



Akiu ^C{GkiuKlJ-^ Guu) - • Gkiu) 

Bklab =C(Gsa,fcA(^, GsaM)) ~ M ' H^^^ Gsa,k) + 

C{GsbdA{p, Gsb,i)) - M • A(m, Gsbd). 

Finally, we find the optimum I for maximizing L(/i,Xi) 
subject to ( [TtI i. ( [Tol l and (fTTT i. This problem is equivalent to 
solving 



max 

I,{tfci|Vfe 



^ ^ {tkluAklu + tklabBklab) 
k.l u,a,b 

s.t. ^tfc, -l,Vfc, ^tfc, = 1,V/, (24) 
( fe 

ife-i = ^ ife-iM + ^ tkiab,y k, I. 

u a.b 

tklu > 0, tklab >0,yk,l,u,a, b. 



Note that the inequahty 

{tkluAklu + tklab Bklab) < tklGkl (25) 



E 

u.a.b 



holds where = max{max„ maxa,6 Sfcjab}- Let us 
call Akiu as the metric for tkiu and i?fe;a6 as the metric 
for tklab- This inequality is tightened when all entries of 
{tkiu,tkiab\y u, a, b} are assigned to zero, except that the one 
with the metric equal to Cki is assigned to tki ■ Therefore, after 
the problem 

r , tkiGki 

{tfc,|Vfe-,i} 

s.t. ^tfe,-l,Vfc, ^tfci = l,V/, (26) 

I k 

tkl > o,yk,i, 

is solved for its optimum solution {i^;|V k, /}, an optimum I 
for (|24] | can be constructed by assigning for every combination 
of k and /, all entries in {tkiuitkiabl^ u,a,b} C I to zero, 
except for the one with the metric equal to Cki to f^^. 

Most interestingly, (l26T l is a standard assignment problem, 
hence {i^;|V k,l} can be found efficiently by the Hungarian 
algorithm, and every entry in {t^;|V k,l} is either or 
1 ITSl . After knowing {i^JV k,l}, the optimum I can be 
constructed according to the way mentioned earlier Finally, 
the corresponding Xi = {I, Pj} is assigned to X^j. Note that 
the Hungarian algorithm to solve (l26t has a complexity of 
0{K^) dgj, meaning that the complexity of finding X^ is 
0(A-3). 

To find /i*, an iterative method which updates ^ with /i = 
[/X— (5(Pt — P(X^))]+ can be used, where (5 > is a prescribed 
step size ifTTll . However, this method converges very slowly, 
since S has to be very small to guarantee convergence. To 
develop a faster algorithm, we first prove that P(X^) is a 
decreasing function of fj, > 0. To this end, suppose fii > 
fi2 > 0. Since Pt — P(Kf^^) is a subgradient of at /i, the 



inequalities d{ni) > d{p2) + /^2)(Pt - ^(Xpa)) and 

d{tJ-2) > rf(Aii)+_(M2 -Mi)(^t -P(X^J) follow. As a result, 
(Mi-/^2)(Pt-P(X^J) ^d(Mi)-rf(M2) > (/ii-M2)(Pt- 
P(X^J) holds, and thus P(X^J) < P(X^J, meaning that 
P(X^) is a decreasing function of > 0. Based on the above 
property, a bisection method can be used to find the n > 
satisfying P(X^) = Pt as /.j*. 

Algorithm 1 The algorithm to solve ( fTSl l. 
compute Gkiu, V k,l,u. 
Mmin ^0; ^inax = 1; computc P(X^„,_^^); 
while P(X^,„^J > Pt do _ 

/^max = 2^max; compute P(Xp„^^); 
end while 

while ^max - Mmin > dO 

if Pt - e < P(X,,) < Pt then 

go to line 12; 
else if P(X^) > Pt then 

Mmin — 

else 

end if 
end while 

compute the S corresponding to X^ and output it as an 
(at least approximately) optimum for (PI). 



The overall procedure to solve (IT&I is shown in Algorithm 
[T] where e > is a small prescribed tolerance. It can be 
shown in a similar way as in |[l5l that the finally produced 
X^ is either equal to (if P(X^) ~ Pt is satisfied), or a close 
approximation (if Pt — e < P(Xp) < Pt is satisfied) for 
X*. Moreover, the indicator variables in X^ are either or 
1. Therefore, the S corresponding to X^ is either optimum or 
approximately optimum for (flST l. After finding this optimum 
S, the optimum subcarrier assignment and source/relay power 
allocation can computed accordingly. It can readily be shown 
that Algorithm [T] has a polynomial complexity with respect to 
K and U. 

IV. Numerical experiments 

Consider the relay-aided downlink OFDMA system illus- 
trated in Fig. 12] t/ = 5 users are served and the users are 
randomly and uniformly distributed in a circular region of 
radius 50 m. To evaluate the maximum sum rate of the adopted 
protocol for each fixed set of system parameters, 500 random 
realizations of channels are generated. For each realization, 
the user coordinates are first randomly generated, and then 
the channels are generated in the same way as in |fT9l . 

When K = 32 and Pt/a^ increases from 15 to 25 dB, the 
average P„ov and Pbcn which are the optimum sum rates for 
the two protocols, respectively, are shown in Fig. |3] It can be 
seen that the average Pnov is always greater than the average 
Pbon- This illustrates the benefit of using the novel protocol. 

When Ptjcr' = 20 dB and K increases from 4 to 64, the 
average are shown in Fig. |4] It can be seen that the 
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Fig. 2. The relay-aided downlink OFDMA system considered in numerical 
experiments. 
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Fig. 3. The average Rnov and -Rbon when K = 32 and PtjcP" increases 
from 15 to 25 dB. 



average is always greater than 1, and the benefit of using 
the novel protocol increases as K increases. This indicates 
that the proposed protocol can better exploit the frequency - 
selective fading than the benchmark one for optimizing the 
subcarrier assignment to maximize the sum rate, especially 
when a big number of subcarriers is used. 

V. Conclusion 

We have proposed a novel subcarrier-pair based opportunis- 
tic DF relaying protocol for downlink OFDMA transmission. 
Note that the proposed protocol truly improves the DF relaying 
itself. It is shown that the novel protocol leads to a maximum 
sum rate greater than or equal to that for the benchmark one. A 
polynomial-complexity RA algorithm has been developed for 
each protocol to maximize the sum rate of all users. Numerical 
experiments have illustrated that the novel protocol can lead 
to a much greater sum rate than the benchmark one. 
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